The fate of orbitons coupled to phonons 
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The key feature of an orbital wave or orbiton is a significant dispersion, which arises from ex- 
change interactions between orbitals on distinct sites. We study the effect of a coupling between 
orbitons and phonons in one dimension using continuous unitary transformations (CUTs). Already 
for intermediate values of the coupling, the orbiton band width is strongly reduced and the spectral 
density is dominated by an orbiton-phonon continuum. However, we find sharp features within the 
continuum and an orbiton-phonon anti-bound state above. Both show a significant dispersion and 
should be observable experimentally. 

PACS numbers: 71.45.Gm,71.38.-k,71.28.+d,71.20.Be 
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I. INTRODUCTION 

In correlated electron systems, orbital degeneracy 
is discussed as an interesting source for very rich 
physicsi^i^. The orbital degeneracy may be lifted 
by a coupling to phonons^ and/or by superexchange 
processes^. The latter give rise to an intimate connection 
between spin and orbital degrees of freedom. Spin-orbital 
models predict interesting ground states such as orbital 
order or quantum-disordered orbital- liquid states, and 
one may expect novel elementary excitations: dispersive 
low-energy orbital waves termed orbitons^. Experimen- 
tally, the observation of orbitons has been claimed on the 
basis of Raman data of the orbitally ordered compounds 
LaMnOgi and RVO3 (R=La, Nd, Y^. In LaMnOg, the 
bosonic orbitons are treated similar to magnons in a long- 
range spin-ordered state, whereas the orbital excitations 
of RVO3 are discussed in terms of one-dimensional (ID) 
fermions equivalent to the spinous of a ID Heisenberg 
chain^i^. 

The role attributed to phonons varies largely. For the 
case of LaMnOs, it has been argued that the orbiton- 
phonon coupling constant g is small, only giving rise to a 
small shift of the orbiton band^-^. The dynamical screen- 
ing of orbitons by phonons and the mixed orbiton-phonon 
character of the true eigen modes has been studied by van 
den Brink^^ using self-consistent second order perturba- 
tion theory (SOPT) in the orbiton-phonon coupling g. 
He interpreted the Raman peaks observed at about 160 
meV in LaMnOa^ in terms of phonon satellites of the 
orbiton band, assuming weak coupling g. In contrast, 
Allen and Perebeinos-'-l argued that the coupling g is so 
strong that the eigen modes can be described in terms 
of local crystal-field excitations ('vibrons'). Then, the 
spectral density consists of a series of phonon sidebands 
(Franck-Condon effect) with a center frequency > 1 eV^^. 
Experimentally, the orbiton interpretation of the Raman 
features observed at about 160 meV in LaMnOs has been 
strongly questioned based on the comparison with in- 
frared dat a^^i^^ , which clearly indicate that these features 
should be interpreted in terms of multiphonons. 



The 'non-local' collective character of the orbital exci- 
tations is relevant if the energy scale of the superexchange 
is larger than the coupling to the phonons. At present, a 
systematically controlled quantitative description of the 
gradual transition from well-defined dispersive orbitons 
at 5=0 to predominantly 'local' crystal-field excitations 
for strong coupling is still lacking. The SOPT treatment 
is valid at small g. Qualitatively, it shows (i) a polaronic 
reduction of the orbiton band width^*'-. due to the dress- 
ing of the orbiton by a phonon cloud yielding a heavy 
quasi-particle, and (ii) a transfer of spectral weight from 
the orbiton band to phonon side-bands, i.e., to a broad 
and featureless orbiton-phonon continuum (see below). 
This suggests that signatures of collective behavior are 
rapidly washed out with increasing g. 

We perform a well-controlled calculation of the spec- 
tral line shape of the orbitons at larger values of g which 
elucidates the collective character of the excitations. Our 
study is based on continuous unitary transformations 
(CUTs) ^"^d^d6 realized in a self-similar manner— in real 
space^". For clarity and simplicity, we restrict ourselves 
to the minimal model in one spatial dimension (ID). 
We find that well-defined, dispersive features appear in 
the orbiton-plus-one-phonon continuum for intermediate 
values of g. Additionally, a sharp orbiton-phonon anti- 
bound state (ABS) is formed above the continuum. Both 
phenomena should allow to observe dispersive signatures 
of orbitons in experiment. 



II. MODEL 

We study the following Hamiltonian in one dimension 

i i 

+ E + E clc. (&I + k) (1) 

i i 

where c and h are bosonic operators that represent 
orbitons and phonons, respectively, w™^, and lo'^^-^ = 
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denote the contributions to the local orbiton en- 
ergy from superexchange and from a static Jahn- Teller 
deformation^^, J the superexchange coupling constant, 
tjph the phonon energy, and g the coupling between or- 
bitons and phonons. The Jahn- Teller energy = 
4(7^/a;ph results from the static distortion of the local en- 
vironment of the transition-metal ion with orbital degen- 
eracy. This distortion lifts the degeneracy such that the 
orbitals are split into a low-lying orbital termed ground 
state orbital and a higher-lying one. The orbital excita- 
tion consists in lifting an electron from the ground state 
orbital into the higher-lying one. This is the effect of the 
creation operator . The superexchange J acts as a hop- 
ping amplitude of the orbiton, giving rise to a dispersion. 

The Hamiltonian ^ is the bosonic version of the well- 
investigated ID Holstein model of coupled electrons and 
phonons yielding polarons, see e.g. Refs. IT9l|20l At T = 
0, Eq. ^ can be studied with a single orbiton so that 
its statistics does not matter and the fermionic and the 
bosonic model are equivalent. 

The dynamic orbiton-phonon interaction (last term in 
Eq. ([l])) corresponds to the creation or annihilation of a 
local distortive phonon if an orbiton is present. This is 
qualitatively the most important term linking the orbiton 
and the phonon degrees of freedom because it represents 
an interaction of one orbiton with one phonon. Of course, 
a hybridization term linear in both the operators of the 
orbiton and of the phonon will also be present. The hy- 
bridization makes orbital effects visible in the phonon 
channel and vice versa. This fact is important for as- 
sessing the possibility to observe phonons or orbitons by 
certain experimental probes. But such a hybridization 
term does not change the character of the excitations 
qualitatively. One may imagine that the bilinear term 
has been transformed away beforehand by a Bogoliubov 
diagonalization. 

Compared to the complex Hamiltonian studied by van 
den Brink for LaMnOa"'^^, our model is stripped to the 
minimum by allowing only for one local (optical) phonon, 
by being one-dimensional, and by omitting the bilinear 
hybridization term between orbitons and phonons dis- 
cussed above. The locality of the phonon ensures that 
all dispersive effects result from the orbital channel. We 
expect that our results will be generic for higher dimen- 
sions also because we do not focus on particular one- 
dimensional aspects. For example we consider a sin- 
gle orbital excitation, not a dense liquid of excitations 
prone to display specific one-dimensional physics such as 
Luttingcr- liquid behavior. For the same reason, we do 
not need to consider the formation of 'biorbitons', the 
equivalent of the bipolarons studied in the context of the 
Holstein model. Hence we are convinced that our model 
contains the generic features, while it is clear and simple 
enough to allow for the controlled computation of line 
shapes. 

In the literature on the Holstein model, the crossover 
from local polarons to large polarons is treated mainly by 
numerical technique s^^'^*^ . The focus has been on tjph < 
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FIG. 1: Dispersion uj{k) of the dressed orbitons for ojph = W 
using the truncation i^V ,e = 9). Inset: Effective band width 
Wc£F as obtained by CUT and by SOFT. For small values of 
g one has 5Wc« oc —g^ for W < Wph and SWeB oc — ff*^^ for 
W = ojph in ID. 
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FIG. 2: Effect of the different truncation schemes on the or- 
biton dispersion aj(fc) for tjph = W . Upper panel: the differ- 
ence SLj{k) = io{^r, e — 9) — uj(^V, e = 9) at constant extension 
e = 9. Lower panel: Suj{k) = Lj(^V,e = 7) - aj(^V, e = 9) at 
constant structure S = ^V. 



W where W is the bare electron (here: orbiton) band 
width (in ID W=J), whereas Wph > IF is reasonable for 
the description of orbitons in transition-metal oxidea^^. 
Then, the continua of different number of phonons do not 
overlap in energy, allowing for the formation of bound 
states below the scattering states or of anti-bound states 
above them. Our calculation elucidates the regime Wph > 
IF and has the merit to provide an explicit interpretation 
of the features found. Thereby, the understanding of the 
nature of the excitations is enhanced. 
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III. METHOD 

We choose the CUT approacbi^ii^ii^ because it is 
an analytical approach which provides an effective model 
that can be understood more easily. The CUT is defined 

by 



(2) 



to transform H{1 = 0) = ifo from its initial form to an 
effective Hamiltonian H{oo) at I — oo. Here / is a con- 
tinuous running variable which parametrizes the unitary 
transformation C/(0 yielding H{1) = U{l)^HU{l). The 
infinitesimal anti-Hermitian generator of the transforma- 
tion is denoted by rj(l). The properties of the effective 
Hamiltonian depend on the choice of the generator. 

In this work a quasi-particle-conserving CUT is used 
which maps the Hamiltonian Hq to an effective Hamil- 
tonian H{oo) which conserves the number of quasi- 
particleai^ii^. This is suggested from the physics of our 
model. The free orbital wave will be dressed by the in- 
teraction with the phononic degrees of freedom. The ele- 
mentary excitations of the interacting model are orbiton 
excitations with a cloud of phonons renormalizing the 
properties of the bare orbital waves. 

The following choice for the matrix elements of the 
infinitesimal generator rj 



Vi.jil) = sgn {qi - qj) Hi^j{l) 



(3) 



in an eigen basis of iJph = '^ph'^ibjb- is appropriate. 
The qi are the eigen values of i?ph-^- Obviously, the flow 
([2]) stops if either Qi — qj or Hi j = for pairs of i and 
J. Given that there is convergence for / — ^ oo, i.e. iJ(oo) 
exists, this tells us that H{oo) is block diagonal, i.e. it 
conserves the phonon number. 

The convergence is not as easy to see generally. But 
assuming an almost diagonal, non-degenerate Hamilto- 
nian with qi < qj <^ iJ^ ^ < Hjj the leading order in the 
non-diagonal matrix elements fulfills 



diH^j = -~sgn{qi - qj) {His - Hjj)H,^j{l) (4a) 
= (4b) 

implying convergence according to 



iJij cx exp(-|iJi_i - iJjj IO 



(5) 



T he g eneral derivations of convergence are given in Refs. 

Mm 

The commutators required for the flow ^ are com- 
puted using the standard bosonic algebra. The trunca- 
tion of the proliferating terms in the flow equation will 



be discussed below. We consider 
H{£) = ^J„cl+„c, -Ka;ph^6l6, (6) 

i.n i 
2,ni ,712 

-I- J r h 

i,ni ,n2 ,«3 

~^ 5Z ( -'^niin2''l+ni'^i+n2^i+«3^i+«4^i -|- h.C.^ 
i,nj 

The only finite amplitudes at are Jo(^ = 0) = 

^o?b + <^b> Mi = 0) = -J/4, and ro,o(^ = 0) = 2g. 
The generator r]{t) is chosen to be 

m = E (rn„n24+n,C,+„,6l-h.C.) (7) 

4, Til ,r!2 

+ \^f2pn3,n4j I^t u\ 5 _hc 



In the following we denote the considered truncations 
by {S, e) with the structure S G {F, V, ^F, '^V} and the 
spatial extension e. The structure S is defined by the 
maximum number of creation and annihilation opera- 
tors appearing in the flowing Hamiltonian, e.g. S = '^V 
for H{t) given in Eq. [HI The spatial extension e is de- 
fined for a given operator c|+„^ c,+„^ &|+„3 . . . 6, by the 
distance between the leftmost and the rightmost local 
operator, i.e. for e=9 only the exchange amplitudes J„ 
with n G {— 4; — 3; . . . ; 3; 4} are finite. We considered 
e e {1, 3, 5, 7, 9} for the numerical evaluation. We have 
found that e—9 and S — \s sufficient to obtain nu- 
merically stable, accurate results (see below. Figs. [U El 
andO. 



IV. RESULTS 
A. Orbiton Dispersion 

In Fig. [Tl the dispersion of the dressed orbitons is 
shown for Wph — W using the truncation (^U, e = 9). 
With increasing g/W the effective band width W^s is 
strongly reduced (see inset). The orbiton becomes dy- 
namically dressed by a phonon cloud enhancing its effec- 
tive mass. Already for g/W=l/2 we find Wcff/W^=0.35 
in CUT. In contrast, the approximate perturbative result 
Wcfi/W—Q.bQ obtained by SOFT clearly underestimates 
the reduction of the band width. For £i/W^=3/4 we find 
Wcff/M^=0.1 in CUT, a factor of 4 smaller than in SOFT. 
This indicates that for g ^ Q.lbW the excitation may be 
regarded as 'local' for practical purposes, i.e. rather as a 
vibron than as a propagating orbiton. 
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The dependence of the orbiton dispersion on the trun- 
cation level is depicted in Fig. [21 The upper panel shows 
the difference du;{k) = Lu{^T,e^9) - Lo{^V,e^9), illus- 
trating the dependence on the structure S. Clearly, 6u}{k) 
increases with g but remains small {6lu/WcS < 1%) for 
all considered values of g. The truncation error result- 
ing from the finite spatial extension e of the operators is 
even smaller, as shown in the lower panel of Fig. [51 This 
reflects the fact that the physics of the model becomes 
more local with increasing g, hence less extended opera- 
tors are sufficient to describe the orbiton dispersion. This 
is underlined by the observation that the truncation er- 
ror related to the extension e is smaller for g/W —1/2 
than for g/W = 3/8. 

B. Spectral Properties 

For the orbiton line shape, we transform the local or- 
biton creation operator 0^°"^ — cj to an effective one by 
the same CUT as applied to Hq. We study all terms with 
one orbiton operator plus up to four phonon operators^^ 

Oil) = ^rcl + E {Bllclhl + C-<6,J (8) 

r ri,r2 
+ ... 

+ J2 iKr^'-'^^'^AbUblbl + ...)+ h.c. . 

Initially, the only finite value is Ao(^ = 0) = 1- Finally, 
at r = 0, there are five types of contributions (omitting 
spatial indices): c^, c^b\ c'b^b^ c'b^b^b^ and c'b^b^b^b^ . 
For each type we sum the moduli squared of all coeffi- 
cients, yielding the spectral weights given in Fig. [3l 

The sum of all contributions has to equal unity. The 
inset of Fig. [3l shows that the sum of the five consid- 
ered contributions is very close to 1 for g/W < 1/2 so 
that higher-order terms are indeed not important in this 
parameter range. For small values of g/W, the spectral 
weight resides almost entirely in the fundamental orbiton 
band, i.e., the eigen mode is a well-defined orbital wave. 
With increasing g/W spectral weight is transferred to the 
orbiton-phonon continua, which reflects the change of the 
excitation character from an orbital wave to a vibronic 
excitation with increasing g. Already for g/W— 1/2 the 
orbiton-plus-onc-phonon band dominates the spectrum. 
This agrees roughly with the result obtained for the local 
limit of both orbitons and phononsiS. 

Now we turn to momentum (k) and frequency (a;) 
resolved spectral properties. The fc-resolved spectral 
weights of the dressed orbiton and of the orbiton-plus- 
one-phonon continuum are plotted in Fig. 4. The depen- 
dence of these two quantities on the truncation level is 
illustrated in Figs. 5 and 6. Again, the differences be- 
tween different truncations are small for all considered 
values of g. As mentioned above, the physics becomes 
more local with increasing g, explaining the larger sensi- 
tivity on the spatial extension observed for g/W=l/8. 




FIG. 3: (Color online) Spectral weights of the orbiton quasi- 
particle band and of the orbiton-plus-n-phonon continua 
(n=l-4) for ujph/W — 1 using the truncation (^V,e = 9). 
Inset: The sum of the weights of the five considered contri- 
butions obeys the sum rule very well. Lines are guides to the 
eye only. 
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FIG. 4: The fc-resolved spectral weight of the orbiton quasi- 
particle band (top) and of the orbiton-plus-one-phonon side- 
band (bottom) for ujph = W using the truncation {^V , e — 9). 

The corresponding /c-resolved spectral densities are 
shown in Fig. 7. The transfer of spectral weight away 
from the orbiton band is largest at the Brillouin zone 
boundary, see top panel of Fig. 4, where the orbiton 
band is energetically close to the continuum, see Fig. 7a. 
For small values of g/W, the spectral weight 'leaks' into 
the continuum where it appears as a broad hump. For 
g/W^l/8, the results of SOFT and CUT for the line 
shape and the spectral weight still agree well with each 
other, see Fig. 7a. With increasing g/W, the effective 
orbiton band width Wcs is reduced so that the sepa- 
ration between the orbiton band and the orbiton-plus- 
one-phonon continuum increases. At the same time, a 
sharp resonance appears within the continuum, display- 
ing a clear dispersion. 

For g/W=l/2, the orbiton-plus-one-phonon sector 



5 



0.003 
0.002 
0.001 

I 
io -0.001 
-0.002 



-0.01 - 



' I ' I ' I ' I ' I ' I ' I ' I ' I ' 
( r,e=9)-( V,e=9) 

\ 

gAV=l/8 s-^^'- - 

/' gAV=l/4 '^x"'^-- 

gAV=3/8 

gAV=l/2 

I I I I I I I I I I I I I I I I I I I 

(V,e=7)-(V,e=9) 




I . I . I . I . I . I . I . I . I . 

1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 
k[7t] 



FIG. 5: Effect of the truncation on the fc-resolved spectral 
weight of the orbiton quasi-particle band for ujpii/W = 1. 
Upper panel: the difference 51 = J(^r,e = 9) - J(^V,e = 9) 
at constant extension e — 9. Lower panel: SI = I^V, e — 
7) - 7(^1^,6 = 9) at constant structure S = ^V. 
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FIG. 6: Effect of the truncation on the fc-resolved spectral 
weight of the orbiton-plus-one-phonon sideband for uiph/W = 
1. Upper panel: the difference 5/ = /(^F, e = 9) -/(V,e = 9) 
at constant extension e = 9. Lower panel: SI = l(^V,e = 
7) - /(^V,e = 9) at constant structure S = ^V. 




dominates the spectrum, cf. Fig. 3. Its spectral weight 
is almost independent of fc, in contrast to the behavior 
observed for smaller values of g/W, see bottom panel of 
Fig. 4. The relative suppression of the orbiton-plus-one- 
phonon continuum close to the Brillouin zone boundary 
reflects the transfer of spectral weight to the orbiton- 
plus-two-phonon continuum. 



FIG. 7: (Color online) Spectral density of one orbiton 
(red) and the orbiton-plus-one-phonon continuum (black) for 
^pb/W — 1 using the truncation (^V,e — 9). An orbiton- 
phonon anti-bound state (green) is shown where it has gained 
substantial weight. (a) g/W=l/8. (b) g/W=l/A. (c) 
g/W^3/8. (d) g/W^l/2. The 5-functions are broadened 
by r = O.OIW. The dashed lines depict the orbiton-plus-one- 
phonon continuum in SOFT for k = — vr. 
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C. Anti-Bound State 

Interestingly, a very sharp dispersive feature is pushed 
out of the continuum to higher energies. Because it is 
built from states with one orbiton and one phonon and 
because it develops from the scattering states in the con- 
tinuum, but lying at higher energies, this feature is iden- 
tified as an anti-bound state (ABS) of one orbiton and 
one phonon. 

The appearance of an ABS is also found in SOPT, 
though at different energies and with different weights 
(not shown) i2i It occurs for all momenta and arbitrarily 
weak couplings g due to the ID van Hove singularities. 
In contrast, the sharp resonance seen within the orbiton- 
phonon continuum at intermediate couplings is totally 
missed by the SOPT which stays broad and featureless, 
see dashed lines in Fig. 7. 

For 17/1^=1/2 the total dispersion - from the peak 
maxima in the continuum at k=Q to the ABS at k--^ 
- is even larger than the renormalized band width Wcs 
of the orbiton band. This surprising result indicates that 
the dispersion of collective orbital excitations may very 
well be observed experimentally even for intermediate 
values of the phonon coupling g. 

The dispersion of the anti-bound state and of the res- 
onance in the orbiton-phonon continuum comes as a sur- 
prise since we started from entirely local phonons. How 
can an ABS of an orbiton and of an immobile phonon 
propagate ? Since in our model ([1]) a single orbiton can- 
not influence the macroscopic number of phonons we can 
exclude that some hopping amplitude of the phonons is 
induced by g. Hence we are led to the conclusion that 
there must be an important momentum-dependent effec- 
tive interaction between an orbiton and a phonon. A pair 
of orbiton and phonon undergoes correlated hopping. In 
the Hamiltonian Hq in ^ before the CUT, this hopping 
can presumably be understood as a virtual process where 
the orbiton-phonon pair is intermediately de-excited to a 
single orbiton which can hop. Momentum-dependent ma- 
trix elements are certainly also present; but they alone 
cannot explain the dispersion of the ABS. 

In the literature on the Holstein model, the appearance 
of sharp subbands which are qualitatively similar to our 
results has been reported for W significantly larger than 
Wph and large values of the coupling constant . How- 
ever, the identification of the anti-bound state by our 
CUT approach is essential in order to understand why 
the dispersion in higher subbands is larger than in the 



elementary orbiton/polaron band. 



V. SUMMARY 

Our work is intended to provide information on the 
nature of orbital excitations in the presence of substan- 
tial orbiton-phonon coupling in principle. We do not aim 
at any particular compound. Therefore, we have concen- 
trated on a simplified model in one dimension without 
hybridization or phonon dispersion. As indicated in the 
discussion following Eq. (jT]), we expect that our findings 
apply qualitatively also to more specific, extended mod- 
els. In two or three dimensions, the line shape at the band 
edges will change and the anti-bound state will form only 
beyond a certain threshold of the orbiton-phonon inter- 
action. The presence of hybridization will imply that the 
orbiton line shape occurs as a weak feature also in exper- 
imental probes coupling directly to the distortions and 
vice versa. Finally, a finite phonon dispersion will con- 
tribute to the mobility of the collective states. In case 
continua start to overlap, effects of finite life times due to 
decay will be observable. These points summarize what 
we expect for the modifications of our findings in real 
systems. 

We investigated the gradual transition from a prop- 
agating orbital wave to a 'local' vibron with increas- 
ing phonon coupling g. We found already for interme- 
diate couplings a substantial reduction of the orbiton 
band width. The orbiton-phonon continuum is not a 
featureless hump as suggested by second order pertur- 
bation theory (Born approximation), but displays a rel- 
atively sharp, dispersive resonance and an anti-bound 
state. For g/W=l/2, the orbiton-plus-one-phonon sec- 
tor carries more spectral weight than the orbiton itself, 
but it also shows the larger dispersion. Thus the collec- 
tive character is seen rather in the orbiton-phonon sector 
than in the fundamental orbiton sector. This may turn 
out as a major advantage for experiments, since phonons 
and orbiton-plus-phonon features appear in separate fre- 
quency ranges, facilitating the correct assignment. Our 
results demonstrate that signatures of collective orbital 
excitations are not limited to compounds in which the 
phonon coupling g is very small. 

We acknowledge fruitful discussion with H. Fehske and 
G. Khalliulin and the financial support of the DFG in 
SFB 608 in which this project has been started. 
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